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Abstract. We present the first determination of the Galactic polarized emission at 353 GHz by Archeops. The data were 
taken during the Arctic night of February 7, 2002 after the balloon-borne instrument was launched by CNES from the Swedish 
Esrange base near Kiruna. In addition to the 143 GHz and 217 GHz frequency bands dedicated to CMB studies, Archeops had 
one 545 GHz and six 353 GHz bolometers mounted in three polarization sensitive pairs that were used for Galactic foreground 
studies. We present maps of the /, Q, U Stokes parameters over 17 % of the sky and with a 13 arcmin resolution at 353 GHz 
(850 yum). They show a significant Galactic large scale polarized emission coherent on the longitude ranges [100, 120] and 
[180, 200] deg. with a degree of polarization at the level of 4—5 %, in agreement with expectations from starlight polarization 
measurements. Some regions in the Galactic plane (Gem OBI, Cassiopeia) show an even stronger degree of polarization in the 
range 10-20 %. Those findings provide strong evidence for a powerful grain alignment mechanism throughout the interstellar 
medium and a coherent magnetic field coplanar to the Galactic plane. This magnetic field pervades even some dense clouds. 
Extrapolated to high Galactic latitude, these results indicate that interstellar dust polarized emission is the major foreground for 
PLANCK-HFI CMB polarization measurement. 



1. Introduction 

The power spectrum of the temperature anisotropies of the 
Cosmic Microwave Background (CMB) have now been mea- 
sured over most of the relevant angular scales (10 arcmin 
to 90 deg, see a comparison of different experiments in e.g. 
IBenoit et al. 2003al and IBennett et al. 2003> . However, CMB 
polarization is only in its experimental infancy. Theoretical 
predictions are rather tight for the polarization effect coming 
from the last scattering surface. Accurate polarization mea- 
surements are not only useful for breaking some degenera- 
cies between cosmological parameters but also for obtaining 
the gravitationnal wave background. Upper limits on polar- 
ization i |Keating et al. 2001] Ide Oliveira-Costa et al. 2003> are 
now superseded by detections by DASI ( Ko vac et al. 2 002 1 
and WMAP (Kog ut et al. 2003) . New results can be expected 
from BOOMERanG 1 , MAXIPOL 2 and other experiments and 
later from Planck 3 . For high frequency CMB measurements 
the most important foreground is certainly the emission from 
Galactic Interstellar Dust (ISD). Submillimetre and millime- 
tre (hereafter submm) emission intensity of ISD can be in- 
ferred from IRAS and COBE-DIRBE extrapolations {e.g. 
|Schlegel et al. 1998) and has been measured on large scales 
by COBE-FIRAS JReachefa/. 19951 |Boulanger et al. 1996| 
|Lagache et al. 1 998) . On the other hand, nothing is known 
on ISD polarization in emission on scales larger than 10 ar- 
cmin., i.e. those precise scales which are the most relevant for 
CMB studies. It is likely that ISD polarized emission is the ma- 
jor foreground for high frequency CMB polarization measure- 
ments. Ground-based observations of submm ISD polarization 
are concentrated on high angular resolution (arcminute scale) 
of star formation regions. Indirect evidence for large scale po- 
larization come from the polarization of starlight in absorption 
(Fos alba et al. 20021 . Goodman (1996) gives a review of the 
measurements and ambiguities in the interpretation of the back- 
ground starlight polarization. In particular, the visible data are 
biased by low column density lines of sight and do not fairly 
sample more heavily reddened ones. Direct submm measure- 
ments are therefore highly required both for Galactic studies of 
the large scale coherence of the magnetic field and in the field 
of CMB polarization, but are rather challenging as they require 
sensitivities comparable to those of CMB studies. 

Archeops 4 is an experiment designed to obtain a 
large sky coverage in a single balloon flight. First re- 
sults on CMB anisotropies power spectrum are reported in 
(Benoit et al. 2003a Be noit et al. 2003bl . Here, we present the 
first results on ISD polarization measurements with Archeops. 
Its large sky coverage strategy is optimized for finding fairly 
strongly polarized sources without any bias on their location. 
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Section |2 briefly describes the instrument and Sect. [3]the 
ground based calibrations on polarized channels. Section |4] 
presents the specific processing applied to the polarized data. 
More specifically subsect. 14.51 presents the inversion method 
applied to determine the Stokes parameters. Section [5] is ded- 
icated to the main results on local clouds and diffuse regions. 
Section|6]assesses the reliability of the results and Sect0their 
physical interpretation. 

2. Description of the instrument 

Light path 




Fig. 1. a) Scheme of an OMT at 353 GHz used for Archeops. 
The back to back horn is heat sunk to the 10 K stage. The box 
containing the polarizer beam splitter and the mirror is on the 
1.6 K stage. The bolometers on the 100 mK stage and their as- 
sociated horns are not shown. The light enters from the top of 
the drawing into the back to back horn (one horn per pair of 
bolometers). One polarization mode is transmitted through the 
polarizer beam splitter to the first bolometer (A, ||), the second 
mode is reflected to the second bolometer (B, ±). This sys- 
tem ensures that each bolometer of the same pair sees the same 
point of the sky at the same time, b) Orientation of the polariza- 
tion for each bolometer of the three OMTs in the focal plane. 
The arrow represents the scan direction. The thick solid lines 
refer to the || direction, the thick dashed lines refer to the ± 
direction. 

A detailed description of the instrument technical 
and inflight performance is given in ( Be noit et al. 2 002 
IBenoit et al. 2003cb : we here provide only a sum- 
mary description. Archeops is a balloon-borne ex- 
periment with a 1.5 m off-axis Gregorian telescope 
described in ( |Hanany & Marrone 2002) . In partic- 
ular, it satisfies the Mizuguchi-Dragone condition 
( Mizugu chi et al. 1978| |Dragone 1982| l in which there is 
negligible cross polarization at the center of the field of 
view. The cryostat contains a bolometric array of 21 pho- 
tometers operating at frequency bands centered at 143 GHz 
(6 bolometers), 217 GHz (8), 353 GHz (6=3 polarized 
pairs) and 545 GHz (1). The focal plane is maintained at a 
temperature of ~ 100 mK using a 3 He- 4 He dilution cryostat. 
Observations are carried out by rotating the payload at 2 rpm 
producing circular scans at a fixed elevation of ~ 41°. Pointing 
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reconstruction is done a posteriori by using a dedicated 
optical stellar sensor made of a 40 cm optical telescope and 
46 photodiodes. Observations of a single night cover a large 
fraction of the sky as the circular scans drift across the sky due 
to the rotation of the Earth. The experiment was launched on 
February 7, 2002 by the CNES 5 from the Swedish balloon base 
in Esrange, near Kiruna, Sweden, 68°N, 20°E. It reached a 
float altitude of ~ 34 km and landed 21.5 hours later in Siberia 
near Noril'sk, where it was recovered by a Franco-Russian 
team. The night-time scientific observations span 12 hours 
of integration from 15.0 UT to 3.0 UT the next day. The 
polarized channels comprise of three quasi-optical modules, 
which are equivalent to Ortho Mode Transducers (hereafter 
OMT, |B0ifot et q/.T33U| |Hiattopadhy ay & Carlstrom 1999) . 
A pair of conjugated bolometers (see Fig.^ is coupled to the 
sky through a single 10 K back to back horn via the OMT. The 
OMT, which is attached to the 1.6 K stage, is made from a 
single polarizing grid mounted at 45 degrees to the horn axis to 
divide the incoming light into the two orthogonal polarization 
modes. One is transmitted to the first bolometer (A), the other 
is reflected towards the second one (B). At any time, the sum 
of the two bolometer outputs measures the total intensity 
while the difference measures the Q Stokes parameters in the 
OMT eigen basis 6 . The three OMT units are aligned along the 
scan direction and have their (A)-polarization axis oriented at 
60° with respect to each other in order to minimize errors in 
polarization reconstruction (Couch ot et al. 19 99 1 (see Fig.[Q. 

3. Ground-based calibration 

Laboratory measurements were performed in order to calibrate 
the transmission and orientations of the polarizers. The cross 
polarization of a single grid was measured to be less than 1 % 
and is neglected hereafter. 

If t is the intensity transmission rate and if || (±) refers to 
the transmitting (extinguishing) direction of a polarizer, then 
K - (t|| + t ± )/2, k = (t|| - t ± )/2 and q = tJt\\t± are the three 
parameters characterizing a polarizer in Stokes formalism. For 
an ideal polarizer, K = 0.5, k = 0.5, q = 0. If S = (/, Q, U) 
is defined with respect to the observation basis (x, y) and de- 
scribes the polarization state of the radiation propagating along 
— z, and if the || direction of the polarizer makes an angle a 
with x, then the transmitted Stokes vector is S' = MS, with the 
Mueller matrix being 



M 



K 



k cos 2a 



k sin 2a 



k cos 2a K cos 2 2a + q sin 2 2a (K - q) cos 2a sin 2a 
t k sin 2a (K - q) cos 2a sin 2a K sin 2 2a + q cos 2 2a , 



In case of an unpolarized incoming radiation Iq, a photome- 
ter placed behind a polarizer receives KIq. When it is placed 
behind two polarizers that are oriented at angles a and <f, it re- 
ceives (K1K2 +£i&2Cos2acos 2<p + k\k2S\Ti2a sm2ip)lQ. In the 
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case of the 353 GHz bolometers, the OMT polarizer is fixed in 
the focal plane with an angle (p. Rotating a calibration polarizer 
(hereafter CP) in front of it and fitting the measured intensity 
as a function of a gives k\ki cos 2tp and k\ki sin 2tp from which 
ip can be deduced. 

We place a box containing a calibration polarizer which ro- 
tates at 1 .5 rpm above the entrance window of the cryostat. The 
is covered with eccosorb to avoid parasitic reflections. It has 
two apertures. In one aperture we place a 77 K thermal source 
made of polystyrene cup filled with liquid nitrogen. The bot- 
tom of the cup is lined with eccosorb. The other aperture faces 
the cryostat. The black body emission is chopped at 13.4 Hz 
against the ambient temperature to enable lock-in detection. 
This set up allows us to determine the position of the grid polar- 
izers in the OMT (see Fig.Q to within 3 degrees. This source 
of error contributes an uncertainty of less than 5%'mQ and U . 

During the ground based preparation of the flight, we 
placed a matrix of 4 x 4 grids of 50 jjm Cu/Be wires with a 
step of 100 11m in front of the 1 m 2 2 Hz modulated thermal 
source (Benoit et al. 2002 1 placed on a hill at ^ 1 km from the 
telescope. This provided a linearly polarized blackbody source 
for an additional pre-flight polarization calibration. We verified 
that the orientation of the grid polarizers in the OMT agreed 
with the laboratory measurements and found that the beam 
shape in the Q and U states agreed with the / beam shape within 
20%. 

4. Polarization data processing 

4. 1 . Standard processing 

The Stokes parameters reconstruction, as well as the cross- 
calibration described above, only apply to clean data associated 
with an accurate pointing. We here summarize the preparation 
of the data described in more details in (Ben oit et al. 20 03d). 

Pointing reconstruction is performed with about 200 de- 
tected stars per revolution and provides an rms pointing accu- 
racy better than 1.5 arcmin. The polarizers angles determined 
from ground calibrations (see Sect. [3} can then be computed on 
the sky. 

The raw Time Ordered Information (TOI), sampled at ^ 
153 Hz, are decompressed, then filtered to take into account 
the AC biasing scheme coming from the readout electronics. 
Cosmic rays, electronic spikes, artifacts and noisy data are de- 
tected and flagged with an automatic algorithm followed by 
visual inspection. Small areas around strong point sources are 
flagged as well. The flagged data representing less than 1.5% 
of the data are replaced by a constrained realisation of noise 
for subsequent detrending and high-pass filtering. The data are 
corrected with a bolometer model for drifts of the instrument 
response due to changes in the background optical loading 
and in the focal plane temperature. Low frequency drifts due 
to airmass and temperature fluctuations of the various stages 
of the cryostat (0.1 K, 1.6 K, 10 K) are decorrelated using 
housekeeping data (altitude, elevation, temperatures). A spin- 
synchronous atmospheric signal remains and prevents us from 
using the Cosmological Dipole for an accurate calibration. In- 
flight observations of Jupiter lead to the determination of the 
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bolometer time constant (used to deconvolve the data stream) 
as well as the photometric pixel beams (with an error less 
than 10 %). The beams are identical within polarizer pairs and 
moderately elliptical, with a minor and major axis FWHM of 
resp. 10.6 and 13.4 arcmin. We assume that in— flight Q, U 
beams are identical to the intensity beam. 

4.2. Filtering 

We briefly describe here the specific post-processing applied to 
the 353 GHz channels. This processing is not specific to polar- 
ization but rather to Galactic studies for experiments that have 
a scan strategy like Archeops or Planck. 

The major noise component that remains after the pipeline 
(as described above) is some low frequency noise. A brute 
force low pass filter applied on the timeline generates ringing 
on both sides of the Galactic plane. The key issue is thus to 
remove the best low frequency baseline without producing a 
significant ringing. To do this, we first mask the Galaxy us- 
ing a SFD template ( Schlegel et al. 1998 1 and use localized 



slowly varying functions (Benoit et al. 2003c I to interpolate 
the Galactic plane and obtain a first estimation of the baseline. 
This estimation is used to perform a noise constrained realiza- 
tion of the timeline on the masked area. Then, an optimized 
low frequency baseline is calculated using wavelet shrinkage 
techniques ( Macfa s-Perez 20 03 1 which allow to remove high 
frequency noise. This baseline is subtracted from the original 
timeline. 




Fig. 2. Variation of the cross-calibration coefficients a during 
the flight relative to the reference baolometer 353K01. The 
coefficients are mostly constant at the beginning of the flight 
(within error bars), but become noisier after UT=21h30. The 
shaded areas show the +lcr values of the coefficients used in 
the present analysis, and the time interval over which they are 
computed. 



4.3. Cross-calibration method 

Since polarization is obtained from differences of measure- 
ments from detectors at various orientations, it is critical that 
they all be accurately cross-calibrated. Any mismatch in this 
cross-calibration automatically generates intensity leaks into 
the fainter polarization mode. We found that the absolute cali- 
brations obtained on Jupiter and on the Galaxy (Sect. 14.41 are 
not precise enough for polarization measurements: typically, 
in order to detect a 5% polarization on the Galaxy it is nec- 
essary to have a cross-calibration accuracy of better than 2%. 
To achieve higher accuracy cross-calibration we have derived 
a method based on inter-comparing the large signal coming 
from Galactic "latitude profiles" from different bolometers. A 
latitude profile is a tabulation of intensity as a function of lat- 
itude where the data at all longitudes is averaged to produce a 
single intensity value in each 2 degrees latitude bin. Because 
latitude profiles average the intensity from all longitudes they 
have a larger signal to noise ratio compared to two dimensional 
maps of the galaxy. Also, it is plausible to assume that for the 
latitude profiles a spatially uncorrected galactic polarization 
signal averages to zero. This last hypothesis is very important 
and various tests to prove its validity are discussed in Sect. 16. 31 
Let si(b), . . . , s„ bal (b) be the ttboi Galactic profiles (b is the 
latitude bin), measured by riboi bolometers. We make the as- 
sumption that all these profiles are identical up to a calibration 
factor ctj : 

Sj(b) = aj!(b) + nj(b), (1) 



with 1(b) a reference profile and n the noise. Making the as- 
sumption of Gaussian white noise, we minimize the associated 
X 2 with respect to the {aj\j=i, nbol and ~s simultaneously under 
the constraint that a\ — 1. The value of a\ is determined using 
the absolute calibration. We have verified that the cross calibra- 
tion does not depend on which of the {ctj} is the constraining 
parameter. 

Once calculated, these coefficients are used to compute Q 
and U maps (see Sect. I4.5i and to check for the presence of 
a polarized signal. If a residual polarized signal is detected in 
some areas of the sky, we remove these areas before making 
the profiles, and recompute the {otj}j=iji M - We have performed 
simulations showing that with these two steps, the correct rel- 
ative calibration coefficients are recovered with a precision of 
better than 2%. We also checked that the choice of the reference 
bolometer is irrelevant. To perform such simulations we used 
Galactic templates provided by an extrapolation at 353 GHz 
of COBE and IRAS data ( |Schlegel et al. 1998| hereafter SFD) 
and included noise and the polarization properties of the instru- 
ment. 

We calculate the cross-calibration factors using the entire 
data and using only periods of 60 minutes and plot their evo- 
lution during the flight in Fig. |2] The variance in the cross- 
calibration increases starting at about 21hr UT. Around this 
time the scans become more tangent to the Galactic plane. We 
attribute the larger variance to the pattern of the scan and to 
noise induced by the atmosphere. Simulations of the scans and 
a 1// noise model partially reproduce the larger variance; the 
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simulations are limited in their capability to simulate the ac- 
tual noise arising from the atmosphere. For the analysis in this 
paper we keep only data from 15h30 to 21h and Fig.[2]shows 
the value of the cross-calibration factors and their standard de- 
viations that are used for the analysis in this paper. The redun- 
dancy map corresponding to this sky coverage is presented on 
Fig -El The way uncertainties in the cross calibration affect the 
degree of polarization is discussed in Sect. [6] 

4.4. Absolute Calibration 

For absolute calibration we use the so-called FIRAS "Dust 
spectrum Maps" data 7 , which are spectral sky maps (from 
100 fim to 4 mm for each 7° pixel) from which the CMB, in- 
terplanetary dust and interstellar line emission have been sub- 
tracted. Each FIRAS spectrum is fitted with a modified black- 
body emission law: S = T(v/vo'fP(v,Td usr ), where j3 is an 
empirical spectral index. This model is then convolved with 
the Archeops bandpass filter. The Archeops data are smoothed 
to match FIRAS beam. The calibration is then obtained from 
a correlation between FIRAS and Archeops Galactic latitude 
profiles, and has an absolute accuracy of about 6%. This affects 
only the absolute values of /, Q, U and neither the degree of 
polarization nor its orientation. A detailed description of the 
calibration is given in ( Lagache 2003 1 iBenoit et al. 2003ct . In 
order not to include polarization effects in this calibration, we 
calibrated the total intensity of each pair of cross-calibrated 
bolometers against FIRAS. All numerical values throughout 
the paper are given in mKRj. A brightness of 1 iiiKrj is equiva- 
lent to 4.36 MJy.sr -1 using IRAS convention (constant v/ v ) for 
Archeops 353 GHz bandpass filter, and 15.4 mKcMB- 

4.5. Inversion method 

For a given direction of observation n, the associated usual co- 
ordinate vectors (e^, -eg) tangential to the sphere are chosen 
as the reference frame to express Stokes parameters (/, Q, U). 
The sign accounts of Galactic longitudes and latitudes on the 
sphere that are clockwise oriented. Let E be the incident elec- 
tric field, E its amplitude and if/ its angle with respect to in 
the tangential plane, defined in the range [0°, 180°], then 

Q = |E ■ e/ - |E ■ e e | 2 (2) 
U = |E ■ | 2 - |E • ef| 2 , (3) 

where the superscript 45 means that the original coordinate 
vectors have been rotated clockwise by 45 degrees. 

In this subsection, the polarizers are assumed to be cali- 
brated against an unpolarized source (cf. subsec. l4~3*l 14.41 . The 
calibrated polarimeter at an angle a with respect to measures 

m(a) — c E 2 cos 2 (a - ip) + n 

= c (I + Qcosla + U sin 2a) + n, (4) 

7 http : //space . gsf c . nasa . gov/astro/cobe/ 
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where the noise n depends on time, a depends on the bolometer 
and on the pixel, and c is the bolometer calibration coefficient. 
In order to make a n plx map, all samples must be taken into 
account to include noise correlations (in time and from pixel to 
pixel) and equation Q is generalized to : 

M = J\S + N, (5) 

where M is the time ordered vector of the n, x tiboi measures, 
S the (3 n p „)-vector Stokes map of the sky, 9\ the pointing 
matrix encoding the pointing information and polarizer angles 
and N the n, x n\, i noise vector. The^ 2 is given by 

X 1 = (M - JIS) T N-\M - &S) (6) 

and is minimized by the solution 

S = (M 7 N- x Wf x 3? 'N~ l M. (7) 
The covariance matrix is then 

E = GT^AT 1 ^)- 1 . (8) 

Solving the linear system Q is one of the recurrent problems in 
CMB studies since the matrices and vectors are usually large. 
In our case, however, the size of the polarized regions corre- 
spond to temporal frequencies where the noise is essentially 
white (in-scan induced noise), and the level of striping in Q 
and U (cross-scan induced noise) is negligible and therefore 
we can use the following simplification. When the noise is 
not correlated from one measurement to another, N is diag- 
onal and the inversion of large matrices can be avoided. We 
therefore consider each pixel individually, compute the (3,3)- 
matrix Ji 1 ' N~ l 3i and the (3)-vector 3? N~ l M. The system of 
equations thus involves small mathematical objects and the in- 
version time is small. 

5. Results 

Once the data are cleaned (Sect. 14. 1> . they are filtered 
(Sect. 14.21 . calibrated (Sect. 14.41 and combined according to 
Eq. l|5} and inverted with Eq. Q to produce maps of /, Q, 
and U. We choose a pixel size of 27.5 arcmin, corresponding 
to HEALpix ( Gors ki et al. 19 98 1 resolution parameter n s id e = 
128. Pixels that have less than 100 detector samples, which cor- 
respond to 0. 1 1 sec mission integration time and a 1 cr I noise 
level of 143 //Krj are blanked. For display purposes the maps 
are smoothed with a 1 deg beam, and these maps are shown in 
Figs. 1415 loY The noise estimate for /, Q, U is obtained through 
Eq.El 

The dispersion of Q and U at high Galactic latitudes is 
found to be ^ 1.1 times larger than their noise estimates from 
the inversion method. This is due to the fact that the noise is not 
perfectly white. In the analysis that follows we use the mea- 
sured dispersion as a measure of the noise and not the lower 
noise estimated from the inversion method. The instantaneous 
mission / sensitivity is found to be about 48 juKRj.sec 1 ^ 2 . On 
average, the 1 cr noise per pixel of 27 arcmin (« s ;d e = 128) is 
found to be 82 //K RJ in intensity and 105 juK RJ in Q and U. A 
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Number of hits per pixel 




Fig. 3. Total number of detector samples per pixels with the 6 bolometers at 353GHz. The map is centered on the Galactic anti- 
center and grid coordinates are spaced by 20°. The Galactic plane is scanned from about 85 to 120 and 180 to 200°. The covered 
fraction of the sky is 17 %. 



I (mKrj) 




2.Q0e-01 ^^^^^^h ^^^^^h 1.00e + 00 



Fig. 4. Archeops / map at 353 GHz in rnKgy smoothed with a 1 degree Gaussian beam. 



statistically significant polarization signal is detected in various 
locations on the galactic plane. Fig.[7]shows a map of the "nor- 
malized squared polarized intensity" defined as the squared po- 
larized intensity (Q 2 + U 2 ) normalized to its variance (<t 1 q+o 2 ij ). 
Twice this quantity behaves has a x 1 probability distribution 
function with 2 degrees of freedom. A statistically significant 
signal appears in regions where the normalized squared polar- 
ized intensity significantly exceeds a value of unity and several 
such regions are detected along the galactic plane. We now dis- 
cuss the results on isolated regions and diffuse medium. 



5.1. Galactic dense clouds 

Here we focus on connected regions in which the normalized 
polarized intensity exceeds the 2<x level on Fig.0 The Stokes 
parameters for these clouds are determined by averaging the 
pixel values of the unsmoothed maps with weights that are in- 
versely proportional to the variance in the pixels. Table^gives 
the values of /, Q, U, p and 6 for these regions. Because the 
Taurus cloud region has been well studied by various instru- 
ments we give its values separately to facilitate an easy com- 
parison. Since p is not a Gaussian variable, we estimate and 
correct for the bias on its determination using Monte Carlo 
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Fig. 5. Archeops Q map at 353 GHz in mK^y smoothed with a 1 degree Gaussian beam. 
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Fig. 6. Archeops U map at 353 GHz in rnKgy smoothed with a 1 degree Gaussian beam. 



simulations. Error bars on p and 6 at 68 % CL are also de- 
termined using simulations. The simulations include the cross- 
calibration errors discussed in Section l4~3*l 



5.2. Diffuse Galactic regions 

In this section we determine whether there is a coherent level of 
polarization on large regions without defining any cloud bound- 
aries. For that purpose and to enhance the signal to noise ratio, 
we divide the Galaxy into 5 deg wide bands along Galactic 
longitude. For each band we construct three latitude profiles 
consisting of the values of /, Q and U as a function of latitude; 
we use binning of 2 degrees in latitude. These three profiles 



are then used to find a unique polarization vector (p, ff) charac- 
terizing the region corresponding to the profile. We avoid the 
bias in the determination of the polarization vector using sim- 
ulations. An example of a profile is shown in Fig. |8j and the 
results are summarized in Tableland in Fig. [9] Coherent po- 
larization levels of few percents are significantly detected up 
to 5 % at the 3 to 4 cr level for several longitude bands, some 
of which include the clouds already discussed in the previous 
section. Even after masking these clouds, a significant coherent 
polarization remains in the same longitude bands. 
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4.00e + 00 



Fig. 7. Map of the normalized squared polarized intensity {Q~ + U )/(Oq + cry). Twice this quantity is statistically distributed 
like ax 2 with 2 degrees of freedom. The 68, 95.4, 99.7% CL of the mapped quantity correspond to 1.1, 3.1, 5.8 respectively. 
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Fig. 8. Scaling of I, Q, U Galactic latitude profiles e. g. for a Galactic longitude range [105, 110]°. The profiles enable to 
constrain diffuse polarization. 



6. Systematics and cross-checks 

It is the first time that a significant detection of polarization 
from the diffuse submillimetre Galactic dust emission is re- 
ported. Before giving interpretation, we discuss the levels of 
possible systematics that can alter the results. 



6.1. Cross checking different methods of polarization 
determination 

To test our Q and U results, we have employed two additional 
techniques to derive their values. Instead of finding a combined 
solution for I, Q, and U using equation @ we determined 
(Q, U) from differences of cross-calibrated pair of bolometers 
by differencing the time-ordered data. Following equation ® 
one can write 



Am(a) = m(a)/c — m(a + n/2)/c' 

= 2 (2 cos 2a + {/sin 2a') + n' . 



(9) 



where c and c' are the calibration constants of the bolome- 
ters (see Eq. 0}. In this method the information on the total 
intensity / is lost. However, the noise power spectrum of the 
difference Am is much flatter at low frequencies than for in- 
dividual bolometers because the differencing scheme removes 
common spurious unpolarized signals such as the atmosphere 
or common gain drifts. 

We apply the map making algorithm outlined in eqs. (|5][7} 
to the difference Am and re-derive Q and U values. The results 
are consistent within one <x with the results reported in Tables 
[Qand|2]for all clouds and galactic profiles. 
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Fig. 9. Summary figure of diffuse Galactic polarization. The intensity (divided by 10) is represented in black and is taken to be 
the average value for -2 < b < 2 in each longitude band. The thin solid line is the same value in 2° wide bands. The direction of 
polarization for every bin is represented below in green, and the length of the dash is proportionnal to the degree of polarization 
p in %. The horizontal error bars represent the width of the longitude bins, which is 5° except for the edge bins. Values are 
summarized in Tab. 13 



In the second technique we simply bin the difference of the 
TOD of a given cloud in 10° bins of the polarizer angle a. The 
binned signal is fit with a function of the form of Eq. l|9} to 
obtain Q and U values. This method does not depend on the 
map making algorithm which was used both as the main anal- 
ysis technique and for the alternate method described earlier in 
this Section. The results of the binning for the Gemini cloud 
(centered on (l,b) = (194.5,-0.9)) are shown in Fig.[l0| but 
have been carried out for all the clouds for which Table^gives 
results. For the Gemini cloud the fit has a X 2 /ndf of 0.93 and 
gives Q = -0.056 ± 0.006 mK RJ , U = 0.020 ± 0.006 mK RJ 
and = 80.2 ± 3.4 deg. The results show that Eq. is a good 
fit to the binned data and that the Q and U values are in good 
agreement with polarization values deduced from the other two 
techniques described earlier. Figure^|also shows consistency 
between the two sets of three photometers (each with two or- 
thogonal bolometers) and between the data taken at two differ- 
ent times during the flight. 



6.2. Consistency between bolometers 

In order to check the consistency between the three pairs of 
bolometers, we computed Q and U with various combinations 
of only two pairs out of the three available. Figure ^] shows 
that all the results are in good agreement and are consistent 
with the values derived using all three pairs. Moreover, the pho- 
tometric accuracy of the / map can be checked against SFD 
template at 353 GHz. We find a consistency within 10 %. 




50 100 150 

Polarizers angles 



Fig. 10. Fit of difference-timelines (filtered differences of 
bolometer outputs of the same pair) as a function of their as- 
sociated polarizer angle for the Gemini cloud (cloud index 6 
in Table 1). This cloud is observed at two different time in- 
tervals and provides useful consistency check, independent of 
the inversion method. Values are averaged in 10° angle bins. 
Different symbols and colors are for the different pairs Red 
squares, light blue circles and blue triangles are for the first, 
second and third pair resp. Empty symbols refer to the first 
time interval, filled ones refer to the second one. The global 
fit has a x 2 /ndf of 0.93 and gives: Q = -0.056 + 0.006 hiKr,, 
U = 0.020 ± 0.006 mK RJ and G = 80.2 ±3.4 deg, in good 
agreement with values deduced from the maps and mentioned 
in Tab. □ 
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Fig. 11. Q (left), U (right) determined on the Gemini cloud (Z 
193°, b - - 1 .2°) with sets of two pairs of 353 GHz bolometers 
out of the three available. Each point is labeled with the pairs 
used to compute it. The solid line is the value determined on the 
map using the three pairs, the dashed lines are the 1 <j errors on 
this value. 



We can also check the cross-calibration process on the 
three intensity maps that can be deduced from the three pairs 
of bolometers. The histograms shown on Fig.ll2lgive the inten- 
sity difference between two pairs normalized by the expected 
noise. The histograms are Gaussian with a standard deviation 
of 1 and without any outlyers. 

In order to test the reconstruction of the noise characteris- 
tics, we also perform a null test with the data themselves by 
randomizing the angle of each polarizer pair at each sample 
before applying the inversion 0. As maps are computed with 
rather large pixels (« sl y e = 128), this effectively cancels out any 
diffuse sky polarized signal as checked on simulations while 
conserving the noise properties. If one considers the values of 
Q and U derived from these "randomized" maps for the ten 
clouds of Table [0 one can form the x 2 /ndf of the hypothesis 
that they all be zero. We find a compatibility at the 90% CL. 



6.3. Uncertainties on cross-calibration 

The validity of the cross-calibration procedure depends crit- 
ically on the assumption that the regions over which signals 
from different detectors are compared are on average not polar- 
ized. This assumption can fail in various ways and here we test 
for these failures. 

If one region which is highly polarized biases the cross- 
calibration coefficients, then the polarization derived in other 
locations in the galaxy should be correlated with intensity. 
Visual inspection of the maps does not reveal an obvious cor- 
relation. For example, the right part of the intensity map (near 
Cygnus), which is the brightest (Fig. [4}, has no polarization 
counterpart. 

A large scale coherent polarization can induce a system- 
atic error in the cross-calibration if the polarizers crossed the 
Galaxy at constant angles. Because of our scan pattern the po- 
larizers rotate as they cross the galaxy. To understand the ef- 
fect quantitatively we carried out simulations where the entire 
galaxy ( Schle gel et al. 1998) was polarized at the 5% level with 
constant orientation. The reconstructed cross-calibration coef- 
ficients, derived using the assumption that the galaxy does not 
have a large scale coherent polarization, were biased at the 2- 
3% level and modified the polarization at the level Ap ~ 1 %. 

The galaxy is most probably not polarized at a constant 
level and orientation. We derive the cross-calibration errors by 
simulating several test maps of the Galaxy using SFD tem- 
plates with a 5% polarization with random orientations in var- 
ious regions. We perform a first iteration of the cross calibra- 
tion assuming that these simulated galaxy maps are not polar- 
ized, and we then reconstruct the polarization map using the de- 
rived cross-calibration coefficients. We perform a second cross- 
calibration after masking the regions that were found to be po- 
larized after the first iteration with a significance level of larger 
than 2 <x. The l<x uncertainty in the extraction of the cross cal- 
ibration coefficients after this second iteration was 2%. The lcr 
uncertainty in the determination of the Stokes parameters in 
the regions that had 5% polarization was less than 1% in /, and 
about 3-5 % in Q and U, depending on the intensity of the 
chosen polarized region. 



6.4. Filtering effects 

Using the above mentioned simulations (Schlegel et al. 1998 1, 
we observe that time domain filtering removes the large scale 
diffuse emission (which broadly has a cosecant law behaviour). 
It represents usually 10-20 %, sometimes 30 % of the total in- 
tensity along the line of sight. By changing the filtering param- 
eters on the real data timelines (e.g. the mask, the frequency 
cut), we observe a similar effect, principally on / and much less 
on Q and U. We derive the systematic error bars (Tabs.^ an d 
13 from the dispersion of the results obtained with these differ- 
ent filterings. The polarized emission characterized by Q and 
U and its orientation are therefore more accurately determined 
than the degree of polarization. 

6.5. Beam effects and time constants 

Beams are found to be nearly identical between the two 
bolometers of a same pair. Nevertheless, slight beam mismatch 
convolved with the Galactic gradient could generate a spurious 
polarization signal. This effect has been estimated to be of at 
most 10 fiK, below our statistical uncertainties. This is also true 
for uncertainties on time constants which are less than 2 mil- 
lisec (i.e. 6 arcmin). The cross-calibration found using Jupiter 
agrees with the cross-calibration using the Galactic profiles to 
better than 1 .5 cr uncertainty for each bolometer. 

7. Interpretation of the results 

The emission of two cloud complexes appear to be strongly po- 
larized at 353 GHz. One large complex is in Cassiopeia (clouds 
1-5 in Tab. ^ with an area of 33 deg 2 . This area includes the 
supernova remnant CasA, although the center is detected in the 
processing as a point source and is not projected. The other 
complex coincides with the southern part of Gem OBI (cloud 
6 in Tab.[Q. Interestingly, the observed part of the Cygnus com- 
plex is not fount to be significantly polarized. 

The orientation of the polarization that we find using the 
galactic profiles is found to be coherent on large scales and is 
also consistent with that found in clouds. Overall, the orien- 
tation is nearly orthogonal to the Galactic plane. It was long 
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Fig. 12. Histograms of the noise normalized difference between the intensity from two pairs of bolometers on the Galaxy 
(b e [-10, 10]). These histograms are compatible with a normal distribution, which shows a good relative calibration recon- 
struction. 



noted from optical polarization studies that neighbouring stars 
had similar polarization directions, with a degree of polariza- 
tion pv more or less proportional to reddening giving an em- 
pirical relation : p v - 0.03Ay ( Whitte TT996l IGoodman 19 96 1. 
The basic explanation is that a large scale Galactic magnetic 
field induces alignment of elongated ISD grains. Starlight po- 
larization measures the projection of the direction of polar- 
ization on the plane of the sky. However, optical polarization 
measurements sample only rather low reddening lines of sight 
and near infrared polarimetric studies yield ambiguous results 
concerning denser clouds (Whittet 1996 Goodman 1996 1. On 
the other hand, submm polarization is free from opacity ef- 
fects and samples all the ISD material along the line of sight. 
The 353 GHz band is nearly on the Rayleigh-Jeans side 
of dust thermal emission, so grains of various temperatures 
should not have very different contributions in different radi- 
ation fields. If the grains that produce visible extinction are re- 
sponsible for submm emission with the same efficiency, then 
an average polarization of at least 3 % (Stein 1966 1 is ex- 
pected at 353 GHz. As shown in Tab. |2 we find a level 
slightly above this figure and much higher in some clouds 
therefore indicating a very efficient grain alignment mechanism 
( |Hildebrand & Dragovan 1995) . Moreover, starlight extinction 
polarization measurements are predicted to be orthogonal to the 
polarized thermal emission. Catalogs of starlight polarization 
have been gathered (Fosalba et al. 2002 and references therein) 
and show a global orientation parallel to the Galactic plane in 
this longitude range, compatible within 20 to 30° with the ori- 
entation of diffuse medium emission as shown in Tab. |2] and 
Fig. |9] If the magnetic field follows the spiral arms, one can 
also expect, as we measure only its projection onto the plane of 
the sky, that some longitudes should have a reduced apparent 
polarization (see Fig. 5 in iFosalba et al. 20 02 1. The very low 
polarization found on Cygnus is in qualitative agreement with 
this prediction as the spiral arm lies along the line of sight in 
this longitude range. 

A coherence of the orientation of polarization between the 
diffuse medium and denser clouds is generally observed, ex- 
cept for the cloud Gl 13.2-2.7. It seems that the global mag- 
netic field that pervades the Galactic plane also goes deeply 
into some denser clouds and is not tangled by turbulence ef- 
fects. However, the degree of polarization may vary by as much 



as a factor two inside the same cloud complex. This probably 
comes from the local variability of the magnetic field. 

The present observations are complementary to the 
far infrared and millimetre polarimetry as reviewed by 
(Hildebrand 1996 1 because here we probe much more diffuse 
lines of sight. 

Although the instrument sensitivity does not allow to mea- 
sure directly high Galactic latitude dust polarization, we can 
extrapolate our results to these regions, assuming that the co- 
herence of the magnetic field and the properties of the ISD are 
similar to the ones in the Galactic plane. It can then be antici- 
pated that dust polarized emission will be the major foreground 
to CMB polarization studies at the level of 10 % of its intensity, 
as anticipated by ( Prune t et al. 199 8 1. The integration along the 
line of sight of various orientations tends to decrease the over- 
all effect of polarization in the Galactic plane, whereas at high 
latitude, this depolarization effect should be smaller. 



8. Conclusions 

Archeops provides the first large coverage maps of Galactic 
submm emission with 13 arcmin resolution and polarimetric 
capabilities at 353 GHz. We find that the diffuse emission of 
the Galactic plane in the observed longitude range is polarized 
at the 4-5 % level except in the vicinity of the Cygnus region. 
Its orientation is mostly perpendicular to the Galactic plane and 
orthogonal, as expected, to the orientation of starlight polarized 
extinction. Several clouds of few square degrees appear to be 
polarized at more than 10 %. This suggests a powerful grain 
alignment mechanism throughout the interstellar medium. Our 
findings are compatible with models where a strong coherent 
magnetic field coplanar to the Galactic plane follows the spiral 
arms, as observed in external Galaxies. 
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Cloud index 



/ (mK RJ ) (stat) (syst) Q (mK RJ ) (stat) (syst) U (mK RJ ) (stat) (syst) 



1.018 ± 0.015 ± 0.085 -0.057 ± 0.021 ± 0.001 



0.111 ±0.017 ±0.008 



0.562 ±0.014 ±0.090 



-0.115 ±0.019 ±0.001 



0.052 ±0.016 ±0.003 



1.498 ±0.011 ±0.105 -0.004 ± 0.014 ± 0.007 



0.113 ±0.013 ±0.014 



0.419 ±0.021 ±0.111 



0.100 ± 0.028 ± 0.013 -0.018 ± 0.027 ± 0.009 



0.994 ± 0.023 ± 0.075 



-0.125 ±0.029 ±0.015 



0.015 ±0.029 ±0.010 



0.820 ±0.011 ±0.113 -0.135 ±0.014 ±0.010 



0.005 ± 0.015 ± 0.005 



0.698 ± 0.004 ± 0.055 



-0.059 ±0.005 ±0.013 



0.011 ±0.006 ±0.004 



0.409 ±0.010 ±0.039 -0.023 ± 0.014 ± 0.002 -0.010 ± 0.013 ± 0.005 



0.271 ±0.006 ±0.065 



-0.001 ± 0.009 ± 0.002 -0.021 ± 0.008 ± 0.006 



0.473 ± 0.009 ± 0.080 



0.001 ± 0.013 ± 0.003 0.001 ± 0.01 1 ± 0.008 



Cloud index 


/ 


b 


Size (deg 2 ) 


p (%) (stat) 


(syst) 


en 


(stat) 


(syst) 





103.0 


1.8 


5.9 


12 r 1 - 8 


± 1.8 


59 


±4.7 


± 1.0 


1 


105.8 


0.6 


7.3 


22.2^ 


±4.0 


78 


± 3.7 


±0.8 


2 


109.7 
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8.8 


'-0.9 


± 1.5 


46 
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±2.0 


3 


113.2 


-2.7 


2.9 


n+6.5 
-6.7 


±9.7 
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±7.5 


±3.2 


4 


113.6 


-1.2 


2.3 


12.3! 2 ; 8 


±2.6 


87 


±6.7 


±2.8 


5 


115.0 


2.4 


5.9 


16-3!!;? 


±3.5 


89 


±3.2 


± 1.1 


6 


193.0 


0.0 


21.6 


8-5^ 


±2.6 


85 


±2.8 


±3.1 


7 


159.3 


-20.1 


18.5 


5-31:1 


± 1.5 


101 ± 


13.6 


±6.3 


8 


165.6 


-9.0 


50.1 


7 9+2.8 
'■ A -2.S 


±4.1 


133 ± 


11.5 


±3.8 


9 


174.4 


-13.6 


21.0 




< 3.4 


23 ± 16.0 ± 


104.7 



Table 1. Top: Stokes parameters of significantly polarized Galactic clouds (above double line) and of Taurus complex (below 
double line). Bottom: Coordinates of the maximum of intensity, measured area, degree and orientation of polarization of these 
clouds. Systematic error bars are derived from the dispersion of the results with different filtering parameters on the timelines (see 
Sect. l6.4> . Last line correspond to a 95 % CL upper limit. Angles are counted clockwise, being parallel to the Galactic plane. 
A brightness of 1 mKRj is equivalent to 4.36 MJy.sr" 1 using IRAS convention (constant vl v ) for Archeops 353 GHz bandpass 
filter. Absolute calibration error of 6 % is not included. 
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Gal. Long 


Range (°) 


/ (mK 


rj) (stat) (syst) 


Q (mK RJ ) (stat) (syst) 


U (mK 


rj) (stat) (syst) 


85 
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1.34 ± 
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040 
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0.001 ± 
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90 
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0.012 ±0.008 


± 0.002 


0.010 ± 
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95 
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0.085 + 0. 


052 


-0.001 ± 0.008 


± 0.000 


-0.019 ± 


0.007 


+ 0.003 


100 


105 


0.79 ± 


0.081 ±0. 


062 


-0.030 ± 0.008 


± 0.001 


0.037 ± 


0.006 


± 0.005 


105 


110 


1.03 ± 


0.076 ± 0. 


062 


-0.030 ± 0.007 


± 0.003 


0.038 ± 


0.006 


±0.001 


110 


115 


1.21 ± 


0.071 ±0. 


031 


-0.033 ± 0.007 


± 0.000 


0.042 ± 


0.007 


± 0.002 


115 


120 


0.62 ± 


0.069 ± 0. 


029 


-0.009 ± 0.008 


± 0.009 


0.030 ± 


0.009 


+ 0.006 


180 


185 


0.52 ± 


0.068 ± 0. 


048 


-0.027 ± 0.008 


±0.011 


-0.007 ± 


0.008 


+ 0.000 


185 


190 


0.59 ± 


0.061 ± 0. 


055 


-0.028 ± 0.006 


± 0.008 


0.002 ± 


0.007 


± 0.004 


i on 


195 


0.82 ± 


0.041 ± 0. 


055 


-0.044 ± 0.005 


±0.011 


0.010 ± 


O.OUo 


i n nn 1 
+ U.UU1 


1 VJ 


197 


0.43 ± 


0.029 ± 0. 


042 


-0.011 ±0.008 


± 0.005 


0.013 ± 


0.009 


i n f\r\i 




Gal. Long. Range (°) 


p (%) (stat) (syst) 


9 (°) (stat) (syst) 












85 90 




0.2+°;!? ±0.1 


4± 


34 ± 15 












90 95 




0.9+°j ± 0.3 


20 ± 


16 ±6 












95 100 




1.6^ ±0.4 


134 4 


15 ±0 












100 105 




6.0«| ±1.0 


65 


±4±2 












105 110 




4.7+o « ± o 6 


64 


±4± 1 












110 115 




4.3+°;* ±0.2 


64 


±4±0 












115 120 




4.8+Jj ± 1.6 


54 


±8±9 












180 185 




5.2 + | : * ±2.6 


98 


±9±3 












185 190 




4.6+);°, ± 1.9 


88 


±8±4 












190 195 




5 4+ - 6 +17 


84 


±4+1 












195 197 




3.5+^ ± 1.6 


65 4 


17 ±8 







Table 2. Top: Stokes parameters of 5° Galactic longitude bands. The intensity is the average over the longitude band and b 
taken in the range [-2°, 2°]. Q and U are scaled from that intensity using the latitude profile fit. Bottom: Degree and orientation 
of polarization for these bands. Angles are counted clockwise, being parallel to the Galactic plane. Due to incomplete sky 
coverage some longitude bands are not quoted. Systematic error bars are derived from the dispersion of the results with different 
filtering parameters on the timelines (Sect. 16. 4l . 



